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Self consistent renormalization theory of itinerant ferromagnets is used to calculate the Curie 
temperature of clusters down to approximately 100 atoms in size. In these clusters the electrons 
responsible for the magnetic properties are assumed to be (weakly) itinerant. It is shown that the 
Curie temperature can be larger than in the bulk. The effect originates from the phenomenon of 
level repulsion in chaotic quantum systems, which suppresses spin fluctuations. Since the latter 
destroy the magnetic order the resulting Curie temperature increases, contrary to expectations of 
the naive Stoner picture. The calculations are done assuming that the energy levels of the cluster 
are described by the Gaussian Orthogonal Ensemble of random matrix theory. 



I. INTRODUCTION 

Nowadays, nanotechnology is one of the most promis- 
ing research areas. By using atoms and molecules as 
building blocks or by down scaling macroscopic objects, a 
large variety of nano-scale objects can be created. These 
nano-scale objects possess interesting physical properties 
substantially differing from the bulk [1] . It is size quanti- 
zation of energy levels that could cause, sometimes, a dra- 
matic difference of optical, electronic or chemical proper- 
ties compared to the bulk. Applications in material sci- 
ence, health care, etc. are offered for these nano-objects 
0. 

Studying clusters containing tens to tens of thousands 
of atoms is a subject of special interest [3]. Such small 
particles demonstrate unusual mechanical, optical, and 
magnetic properties and chemical activity. The latter is 
crucial, for example, for their use in catalysis processes. 
However, the precise mechanisms behind the physical and 
chemical properties of the clusters are not yet fully un- 
derstood pp. For new applications and improvements 
of existing ones it is extremely important to understand 
these physical principles. 

The study of magnetic properties is one of the methods 
to gather information about the cluster's behavior. 
Namely, the magnetic properties are sensitive to atomic 
and electronic structure, quantum size effects, surface 
to volume ratio, and symmetry. Several experiments 
demonstrate unusual magnetism of metallic clusters. For 
example, it was found that the clusters of 3d transi- 
tion metals iron, cobalt and nickel, show enhanced fer- 
romagnetism in comparison with the bulk [3H5]. Fur- 
ther, rhodium clusters appeared to become ferromagnetic 
while in the bulk rhodium is paramagnet [TJ. For man- 
ganese, a change of magnetic order from fcrro- to anti- 
ferromagnetism occurs during the transition from cluster 
to the bulk regime [8]. 

In general it is a complicated matter to understand 
the physical principles behind the magnetic properties of 
metallic clusters. Within the picture of localized mag- 
netic moments described by the Heisenberg model one 
could expect that the smaller the cluster size and, thus, 



the larger the ratio surface to bulk, the weaker the mag- 
netism is, due to a lower coordination of surface atoms. 
For an itinerant electron system the influence of down- 
scaling is not immediately clear. For certain small cluster 
sizes so called shell effects give rise to an increased den- 
sity of states near the Fermi level [S] and thus, due to 
the Stoner criterion |10j . an increased tendency towards 
ferromagnetism. Besides an influence on the density of 
states, the size could also have an effect on the thermal 
spin fluctuations responsible for the reduction of mag- 
netic order, whose role is crucial for the understanding 
of the magnetic properties at finite temperatures |10j . 

Note that in ultra thin films finite size effects also cause 
interesting physical features such as changes in the Curie 
temperature [TTJ [T^] . For example in Gd ultra thin films 
an increased surface Curie temperature is observed com- 
pared to the bulk [TT] . 

The theory of itinerant electron magnetism is, in gen- 
eral, very complicated and still controversial (for a re- 
view of some modern approaches, see Refs. 13, 14 and 
references therein) . For a special case of weak itinerant 
ferromagnets close enough to the point of the Stoner in- 
stability a very successful self consistent renormalization 
(SCR) theory has been developed [TUHH]. This theory 
clarifies the role of spin fluctuations in finite-temperature 
magnetism and allows us to calculate the Curie tempe- 
rature. Here we use this theory to study the size effects 
on the Curie temperature for metallic clusters, down to 
approximately 100 atoms in size. It appears that suppres- 
sion of the spin fluctuations leads to an increase of the 
Curie temperature when going to smaller cluster sizes. 



II. THE SELF CONSISTENT 
RENORMALIZATION THEORY OF WEAK 
ITINERANT FERROMAGNETISM 

For a convenience of readers, we present in this sec- 
tion a brief overview of the SCR theory, following Ref. 
ITOl We deal with the simplest model of itinerant elec- 
tron magnetism, that is, the Hubbard model with the 
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Here tji is the transfer integral between the Wannier or- 
bitals at sites j and I, e(fc) is its Fourier transform, that 
is, the band energy dependent of the wave vector fc, au a 
is the annihilation operator for the Bloch state, ci a is 
the annihilation operator for the Wannier state, n$\ is 
the occupation number operator for the Wannier state at 
site j with spin-up, U is the Coulomb interaction energy 
between two electrons on the same site and I = M- with 

iVo 

No the total number of atoms. Thus, the first term of 
the Hubbard model is the hopping term or tight-binding 
model and the second term is the on-site Coulomb inter- 
action term. 

The simplest way to approximate the interaction term 
is by using the Hartree-Fock approximation, which leads 
to the familiar Stoner theory: 

njf n il n 3t ( n ii) + n il ( n 3t) ~ ( n 3t) ( n 3l) ■ ( 2 ) 

Here the (. . .) denotes the statistical average of the quan- 
tity inside the brackets. For the case of a ferromagnet, 
the statistical average of the occupation number is sup- 
posed to be site-independent: (rij-f) = (n-j-) = n-j-. Using 
the approximation of ([2| for Eq. ([!]) and by rewriting 
this expression, leads to 
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Here N is the total number of electrons, M is the mag- 
netization and Eke is the single electron energy given by 



Eka = £fc + crA, 
A = IM + fi B H 



(4) 



with a equal to 1 for spin-up and to —1 for spin down, 
fi b the Bohr magneton and H an external magnetic field. 
Further, 2A is the band splitting between the spin-up and 
spin-down bands. This ferromagnetic state is stable at 
temperature T = if Ip(Ep) > 1 (Stoner criterion) with 
p(Ep) the density of states (per whole system per spin) 
at the Fermi level. 

From the results of Eqs. ([3| and Q it is possible to 
derive the thermodynamic properties of the system in the 
Stoner theory [10] . The Curie temperature, which is of 
main interest here, can be calculated from the divergence 
of the static magnetic susceptibility of the system, which 
within the Stoner theory is given by 
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Here xo is the static magnetic susceptibility of a non- 
interacting system (i.e. with 1 = 0). For the bulk xo is 
given by 



Xo (T) = p(E F ) 
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and the accent corres- 



with R = 

\ f J \f J E=Ef 

ponds to the derivative to E. 

Thus, from Eq. ^ and the divergence of the static 
magnetic susceptibility the Curie temperature Tq follows 
from 

l-7Xo(T c ) = 0. (7) 
For the case of weak itinerant ferromagnets when 

< Ip(Ep) - 1 < 1. (8) 
the Stoner Curie temperature ^ is estimated as 

T$cxE F y/lp(E F )-l. (9) 

It is well known that Curie temperatures predicted from 
this equation are too high compared to the experimental 
data [TU1 [T3] . The reason is a neglect of spin fluctuations 
that allow magnetic excitations at lower energy: it is 
easier to rotate the spins than to change their length like 
in the Stoner theory. 

The SCR theory that takes spin fluctuations into ac- 
count is formulated in the following way [10] (an alter- 
native approach based on diagram technique has been 
developed in Ref. 13). The transverse dynamic magnetic 
susceptibility x~ + (q,u) can be formally represented as 
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Here Xo + { c t, UJ ) is the transverse dynamic magnetic sus- 
ceptibility for non-interacting electrons with spin split 
bands, q is the wave-vector and w the frequency. In gen- 
eral the problem is to find Xmi{q, of Eq. ( 10 1 to make 
the expression for the dynamic susceptibility exact. For 
this purpose the total free energy F of a system is ex- 
pressed in terms of the exact transverse dynamic suscep- 
tibility using the Hcllmann-Feynman theorem: 

F = F + AF 
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Here Fq is the free energy of a non-interacting system. 
Then, one can use the self consistency condition to find 
Xmi(Qi w). This means that the static magnetic suscep- 
tibility calculated from the free energy of Eq. (11) (via 



J^-p = - ) is equal to the static long wave-length limit of 
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the exact transverse dynamic magnetic susceptibility of 
Eq. ((TO ). 



Of course, this condition is not enough to find the 
whole function A. However, for the weak itinerant sys- 
tems where the condition ((8l) is fulfilled one can assume 



that A can be considered in Eq. (11) as a number inde- 



pendent on M, I, u> and q; a formal justification has been 
given in Ref. 15. As a result, A is given by the expression 
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Here a = ixo, /a/ = X W' w )/Xo, d is defined as 
Xo/x = ad=l-a+ X(T,d) and X o = Xo + (0,0). At 
the Curie temperature d = and M — 0, which makes it 
natural to expand /m in terms of small u) and q (u <C -Bp 
and g <C &f respectively with fcjr the Fermi wave- vector) 
and to approximate the nominators of Eq. ( 12 ) by their 
static and long wavelength limit. Thus, at the Curie 
temperature /m can be approximated for the bulk by 



/o = l- Aq 2 + iC- 



(13) 



where A and C are constants (depending on the shape 
of the Fermi surface), and the subscript zero refers to 
M = [TU]. 

By using the condition of divergence of the static mag- 
netic susceptibility, the Curie temperature for a weak 
itinerant system is given by the equation 



l-/Xo(Tc) + A(T c ,0) = 0. 



(14) 



Here, it can be noticed that it is actually A which takes 
the influence of the spin density fluctuations on the Curie 
temperature into account. For further calculations it is 
convenient to separate A in the following two parts: 



A(T,0) = A + A 1 , 
1 f°° 

X o = 7T du;sign{w}Im{G(u;,0)} , 



(15) 



Ai 



1 f°° 2 
— / rfcjsignjo;} Im{G(^,0)}. 

Z7r J -00 e t 1 



The function Ao is the temperature independent part, 
which merely gives rise to a shift of the Stoner criterion 
at T = 0. It can be simply considered as a renormal- 
ization of the Stoner parameter /. Strictly speaking, Ao 
does depend on temperature via Xo> but its temperature 



dependence can be neglected compared to that of Ai. 
Namely, one can show that Ai cx T 4 / 3 compared to the 
T 2 dependence of xo PU] • Important is that due to this 
temperature dependence of Ai the Curie temperature is 
effectively lowered compared to Stoner theory: 



T c oc E F [Ip(E F ) - 1] 
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In other words, the dynamics of the spin density fluc- 
tuations is crucial for the correct description of the Curie 
temperature. 



III. SIZE DEPENDENT ENERGY LEVEL 
DISTRIBUTION 



The SCR theory described above will be used further 
to calculate the Curie temperature of metallic clusters 
as a function of their size. Note that the transverse dy- 
namic magnetic susceptibility of a non- interacting system 
is in fact the only input required for this calculation. For 
clusters this quantity substantially differs compared to 
the bulk due to the differing energy level spectrum. The 
exact calculation of the energy spectrum of a cluster with 
a given size is actually impossible due to certain random- 
ness of its shape. To overcome this problem, the random 
matrix theory [16 20 is typically used. 

In principle it is possible to calculate the well defined 
energy levels of an individual small particle. For the case 
of a perfect metallic sphere this is actually quite obvious. 
One obtains an energy level spectrum consisting of highly 
degenerate energy levels, where the separation between 
the groups of the energy levels is proportional to with 
L the diameter of the sphere and the high degeneracy is 
due to the high geometric symmetry. It is worthwhile to 
stress, however, that although the nonzero level splitting 
is proportional to -^2 , the average energy level separation 
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around the Fermi level follows the well known ^ oc 
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L 3 



proportionality (with V the volume). 

Obviously the situation of the perfect sphere is very 
special. However, there are examples of clusters with 
high geometric symmetry [5]. These are quasi-spherical 
clusters with so called "magic" numbers N = N m , where 
N is the number of atoms. For example N m — 13, 55 
and 147 for the Mackay icosahedron, where each "magic" 
number corresponds to the right number of atoms so that 
a spherically shaped cluster can be formed by packing the 
Mackay icosahedrons in the proper way. Characteristic 
for the energy level spectrum of these "magic" clusters is 
the shell structure. These are highly degenerate or close 
groups of energy levels, which causes the energy level 
separation near the Fermi level to be rather small. The 
effects of this shell structure on the electronic pairing 
in superconductors has been discussed in Ref. [SI For 
the Curie temperature of the magic clusters the smaller 
energy level separation around the Fermi level could have 
important consequences, too. 

Here a generic case will be considered only, mean- 
ing that the situation of the highly geometrically sym- 
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metric clusters is excluded. It is assumed that there 
are uncontrollable atomic surface irregularities, which 
are sufficient to split apart this large degeneracy of 
the energy levels. Further, it is assumed that the 
clusters are large enough to satisfy the proportionality 
8 = 1/ p(Ep) oc 1/V, where 8 is the average energy level 
spacing around the Fermi level and V is the volume of the 
cluster |17H19j . Thus, considering an ensemble of clus- 
ters of the same size, they will differ in their energy level 
spectrum due to the uncontrollable surface, but have the 
same average energy level spacing around the Fermi level. 

Gor'kov and Eliashberg (2JJ were the first who recog- 
nized that this situation is similar to the interpretation 
of nuclear energy level spectra discussed by Wigner and 
Dyson |20j . The idea was to circumvent the unknown 
and complex interactions between the nucleons by us- 
ing a statistical description leading to an energy level 
distribution. To be more specific, it was assumed that 
the eigenvalues of a random matrix could be taken as a 
model for the energy levels of a complex nuclear system. 
This means that an ensemble of possible nuclear systems 
corresponds to an ensemble of random matrices. Impor- 
tant to remark is that the randomness of each matrix 
is restricted, because they must possess certain transfor- 
mation properties imposed by the symmetries that each 
individual Hamiltonian is supposed to have in common. 
Then, depending on the imposed symmetry properties 
different energy level distributions can be derived |20j . 

In the same manner as described above the uncontrol- 
lable surface irregularities can be interpreted as giving 
rise to a random matrix treatment. For metallic clusters 
the transformation properties of the random matrix are 
determined by the magnitude of the spin-orbit coupling 
H so and the magnitude of the Zeeman splitting in an ex- 
ternal magnetic field 2fisH compared to 8. For example 
for small H so and small 2/ibH compared to <5 the ma- 
trix must have respectively rotational and time-reversal 
invariance (T6ll20j . This example corresponds to the so 
called Orthogonal ensemble, which is used in this work. 
The other possible ensembles are given in Table [I] 



for the cluster system [151 HH HO] . The result is 
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(17) 

where Pm is the probability to find a certain energy level 
spectrum, 7 = 1,3 and 4 correspond, respectively, to 
orthogonal, unitary and sympletic ensemble, Cjj follows 
from the normalization condition and k is an ensemble 
dependent constant |16j . From the product in equation 
(171 the level repulsion can be clearly seen. 
A very important quantity is i?(|_E|), which can be 



derived from the energy level distribution (17 1 and gives 



the probability to find two energy levels separated by 
an energy E independent of the number of energy levels 
in between them [T§]. For the orthogonal ensemble this 
distribution is given by 



R{\E\)=l- 




sin(7ra:) 



(18) 



dx. 



This expression will be used in the following Section for 
the calculation of the transverse dynamic magnetic sus- 
ceptibility for a cluster. 

Until this point no comments were made about the cor- 
rectness of the random matrix theory or the assumption 
that the energy level distribution in an irregular cluster 
is universal and depends only on the symmetry class. At 
this point one can say that it is still a hypotheses that 
needs to be tested more in order to reach complete un- 
derstanding of the situation. However, at the moment 
there are many experiments that appear to confirm this 
theory. Q7H3I] 



IV. RESULTS AND DISCUSSIONS 



TABLE I: overview of different Hamiltonian symmetries rele- 
vant for energy level distribution 



Probability distribution 


Magnetic field 


Spin-orbit coupling 


Poisson 


Large 


Small 


Orthogonal 


Small 


Small 




Small 


Large (even particles) 


Unitary 


Large 


Large 


Sympletic 


Small 


Large (odd particles) 



The Poisson ensemble is typical for systems with a reg- 
ular classical motion, there is no level repulsion in this 
case. In the case of chaotic system (three other ensem- 
bles in Table [I]) the probability to find two levels with 
close energies is suppressed. By taking a proper average 
over the ensemble of the random Hamiltonian matrices, 
one can obtain an energy level distribution satisfactorily 



In this paragraph the Curie temperature of clusters as 
function of their size is calculated. From Eq. (14) it is 



clear that this size dependence could come from xo an d 
A. First the size effect on xo will be considered, second 
that on A and finally the resulting effect on the Curie 
temperature. 



With the use of the probability distribution (17 1 it is 
possible to calculate the static magnetic susceptibility 
for a cluster system. In Ref. 21 an interpolation scheme 
between the regimes T -C 8 and T ^> 8 (or bulk) has been 
suggested for which both well developed approximations 
exist [IS] [22] • The important result of these calculations 
[17l [TBI [22] , which will be used later on, is that already for 
T > 8 (we do not mean strong inequality here!) the static 
magnetic susceptibility of a cluster can be approximated 
by that of the bulk. 
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Before A can be calculated, an expression for the trans- 
verse dynamic magnetic susceptibility of a cluster system 
has to be found. We will follow analysis originally pro- 
posed in Ref. [22for the case of optical polarizability (fur- 
ther this result has been slightly corrected, see for review 
Ref. ITTj) . The starting point is the general expression for 
the transverse dynamic magnetic susceptibility, 



X 
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uj 



iO 



(19) 

where \v) and \ are eigenstates of the system and n(E) 



is the Fermi function. Equation ( 19 ) accounts for a single 
particle. Therefore, for the Orthogonal ensemble under 
consideration, one has to average over the energy level 
distribution given in equation (|17[) (with 7 = 1). One 



can then derive that this comes down to multiplying each 



term in the sum of equation ( 19 ) by 



R(\E V -E^\) 
8 2 



dE,,dE„ 



For later convenience the expression for the transverse 
dynamic magnetic susceptibility of an Orthogonal ensem- 
ble of clusters Xo + ( UJ > l) can De written as 



dE fl dE u S(E - E^5(E' - E v ) 
2 n(E')-n(E) R{\E' - E\) 



Xo + (u,q)= J J dEdE' 

x \(v\e^- F \n) 



E -E> -uj + iO 
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For uj <C Ep, which will occur naturally for the cal- 
culation of A at the Curie temperature (Section II), it 
was shown in Ref. [21] that the matrix element is approxi- 
mately energy independent leading to a separation of the 
q and uj dependencies: 



dEdE 



, n{E')-n(E) R(\E' - E\) 



E — E' 



iO 



S 2 



dE^dE v 5{E F - E„)S(E F - E v ) \{u \e^ ? \ y)\ 



(21) 



An accurate computation of this g-dependence or 
matrix element is difficult. However, for q c <s^q-^kp 
(q c oc 1 /R is the inverse average radius of the cluster) the 
g-dependence can be approximated by that of the bulk, 
because classical trajectories of electrons in this case are 
mainly like in the bulk, with rare reflections at the bor- 
der. Thus, within this q-regime the real and imaginary 
part of the dynamic susceptibility are in highest order 
proportional to, respectively q 2 and 1/q (Eq. (13)). [TO] 

For the w-dependent part of equation (21) it can be 
shown that it equals the complex function 



important consequences, one for the evaluation of Eq. 
( [12] ) and the other for Eq. ( [22] ) . Here the former will be 
discussed first, because the latter will follow from it. 



In this regime (q c <^.q<^kp) the sum over q in Eq. (12) 
can be replaced by an integral, because the integrand is a 
smooth function with a maximum at approximately q oc 
uj 1 / 3 , and relevant frequencies are of order of temperature 
(here we use units H = fcg = 1). For the remaining 
part of the sum it is assumed that no singularities occur. 
Therefore, it is proportional to 1/R and can be neglected 
for large enough cluster sizes. 

We will thus restrict ourselves only to the case of not 
too small clusters, where uj 3> 5 for relevant uj, other- 
wise, the discreteness of q- vectors for spin fluctuations is 
essential, replacement of sum over q by integral is impos- 
sible, and the problem should be solved numerically for 
a given particular shape of the cluster. In this limit one 
can show that the w-dependence of the real part of Eq. 
( 22 ) can be neglected compared to that of the imaginary 
part. Then, the following expression for the transverse 
dynamic magnetic susceptibility (normalized to the static 
susceptibility) of a cluster system can be obtained 



f (uJ,q) = l-Aq 2 + l C- A2{ ^. (23) 
q z 

Here A and C are the same constants as for the bulk 
(equation (13)), z = 2ttuj/8, the zero in the subscript of 
the function / corresponds to M — (Curie temperature 
is found from divergence of susceptibility in paramagnetic 
regime) and A2 is the imaginary part of the function A(uj) 
( 22 ) given by 



A a {z) =2t? 



2 sin 2 (77) 
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V 

" sin(i) 
t 



dt- 



where r\ — z/2. At z— >oo A^iz^jz 
tion to the bulk (Figure [lj . 



d 1 sm(?7j \ 
dr) \ r) J ' 
->1 giving extrapola- 




te 2 f R(\E\) 

A(uj) = 2— / — ——rdE. (22) FIG. 1: The function A 2 /z as a function of z. 

J E l — [uj + ivy ' 

Besides the approximation of the matrix element, the At this point it is possible to estimate the cluster sizes 

restriction of the q regime (q c <tiq-^kp) has two other for which the above described approximations are valid. 
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For this purpose our condition oj 5 and u> oc T can be 
used, where the last proportionality can be easily derived 
from the calculation of the temperature dependence of Ai 
(equation (15)). Then, using an estimation 5 ~ Ep/N 
with TV the total number of electrons for the situation 
where Ep = 10 4 K and, say, the value T c = 20K typical 
for weak itinerant ferromagnets JU] results in N > 500. 
Thus, for five d-electrons per atom this would lead to 
the condition that the above described considerations are 
valid for clusters containing approximately more than 100 



atoms. 

As was mentioned, in the case T ^> S the static mag- 
netic susceptibility can be approximated by that of the 
bulk (Eq. ([6])) [22 j. This means that the important size 
dependent contribution to the Curie temperature must 
come from the Ai term only. It can be derived by substi- 
tution of Eqs. (JsJ) and (23 1 into Eq. ( [12] ) and by approx- 
imating the numerators of Eq. (12 1 by their static long 
wave-length limit, Ai (here called Xduster), that is, 



^cluster 

k < 



.1/3 



sin 



'ttxT\ 



TTxT 

s 



sin(t) 



dt 



cos 



sin {^f) 



1/3 



TTxT 

S 



' TTXT ^ 

v 5 I 



(25) 



It is justified since for the case of weak itinerant ferro- 
magnets the region of small q and w is dominant in the 
integral fity . 

The result for X c i us ter is presented as a function of Tj b~ 
in Figure [2j 
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FIG. 3: (color online) The function Xbuik (or 5 = 0) and 
^cluster for S — 0.1 and 1 is plotted as function of temperature. 
Here S is given in units of K and A cr is the critical value of 
A, where the net magnetization becomes zero. 



FIG. 2: (color online) The function A c i„ ster /T 4//3 as a function 
of T/S, indicating the increase of the influence of spin fluc- 
tuations at higher temperatures and particulary larger sizes. 
The question mark indicates the applicability limit of our as- 
sumptions. 

Obviously, for <5— >0 the constant bulk value 
(y%j§- — B) is reached. Further, it is important to 
notice that a cluster system for a fixed size (or 5) has 
a larger temperature dependence compared to the bulk 
which leads to an increase of the Curie temperature. 
For smaller clusters the enhancement of the Curie 
temperature will be larger. 

It is instructive to show the data in a slightly different 
way, where A is plotted as an function of temperature 
for different sizes (5 = 0, 0.1 and 1.0 in units of K) as 
can be seen in Figure [3j To summarize, the results for 
the Curie temperature of the cluster normalized to the 



Curie temperature of the bulk is plotted in Figure [4] as a 
function of the average energy level spacing normalized 
to the Curie temperature of the bulk. The increase of 
the Curie temperature for smaller cluster sizes is caused 
by the increasing suppression of spin density fluctuations 
for smaller sizes. In Section II it was shown that the 
function A takes into account the influence of the spin 
fluctuations on the Curie temperature. 

To verify this theory quantitatively experiments on 
weakly itinerant clusters should be performed. We 
believe, however, that qualitatively our conclusion is cor- 
rect even beyond the limit of formal applicability of the 
SCR theory (equation Q). The level repulsion should 
suppress spectral density of spin fluctuations at small 
frequencies, and this should be the main effect in the 
temperature dependence of magnetic properties, these 
statements being quite general. The experiments on Fe, 
Co and Ni mentioned in the Introduction j3HS] seem to 
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FIG. 4: (color online) Curie temperature of a cluster system 
normalized to the Curie temperature of the bulk as a function 
of the average energy level spacing also normalized to the bulk 
Curie temperature. 

be in agreement with our conclusion. 

It would be interesting to improve this theory further 
for smaller cluster sizes. For this purpose the discrete- 
ness of the energy spectrum should be explicitly taken 



into account and the influence of the static susceptibility 
will become important. For clusters even smaller so that 
8 oc 1/V is not applicable anymore random matrix the- 
ory will fail. In this regime probably ab initio approaches 
are the only way out. 

To conclude, small particles of itinerant magnets show 
an increase of their Curie temperature when reducing 
their size, in a clear contrast to what one would expect 
for a localized picture. Such enhancement of magnetic 
stability originates from the size induced renormaliza- 
tion of electronic states leading to a suppression of spin 
fluctuations, and may open interesting perspectives for 
application of such systems in nanotechnology. 
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